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COLLISIONAL  EMITTANCE  GROWTH  IN  THE  PARTICLE 
SIMULATION  OF  FOCUSED  BEAM  TRANSPORT 

Recent  work  on  space  charge  limited  transport1  2 has  revealed  the  presence  of  instabilities 
in  a beam  with  a Kapchinskij-Vladimirskij3  (K-V)  distribution.  Computer  simulations  have 
shown  that,  in  many  cases,  these  instabilities  saturate  with  relatively  little  emittance  growth4. 
Because  of  the  importance  of  space  charge  limited  transport  to  the  heavy  ion  fusion  program, 
several  simulation  efforts  have  been  undertaken5'6.  Since  it  is  impractical  to  calculate  the  orbits 
of  a number  of  particles  comparable  to  the  actual  number  in  a physical  particle  beam,  it  is  the 
purpose  of  this  report  to  discuss  the  effects  due  to  using  a much  smaller  number  in  the  com- 
puter simulation. 

In  an  ensemble  of  charged  particles  a force  is  exerted  on  each  particle  by  all  the  others  in 
the  ensemble.  It  is  therefore  possible  to  sum  up  the  forces  on  each  particle  due  to  all  the  oth- 
ers and  use  this  force  to  calculate  the  particle  acceleration  and  then  integrate  to  find  its  orbit. 
The  number  of  operations  in  the  approach  increases  as  the  square  of  the  number  of  particles, 
n\ 

On  the  other  hand,  it  is  possible  to  exploit  the  long  range  nature  of  the  coulomb  interac- 
tion to  design  an  algorithm  which  instead  of  scaling  as  n2  scales  as  />.  When  the  number  of  par- 
ticles in  a Debye  sphere  is  much  greater  than  unity,  an  ensemble  of  charged  particles  interacts 
not  through  the  nearest  neighbor  forces  but  through  the  collective  electric  fields.  Numerical 
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Poisson’s  equation  to  find  the  corresponding  electric  field  on  this  grid.  These  grid  fields  are 
then  used  to  move  the  particles. 

A previous  report'  dealt  with  the  effects  of  varying  various  numerical  parameters  in  this 
approach.  Its  emphasis  was  however,  on  the  effects  after  a short  time  and  the  number  of  parti- 
cles was  chosen  large  enough  so  that  only  minimal  effects  were  seen  by  factor  of  two  changes  in 
their  number.  Since  the  binary  interaction  approach,  though  impractical  for  large  number  of 
particles,  can  be  more  efficient  for  a small  number  of  particles,  it  is  important  to  determine  how 
many  particles  are  actually  needed  and  perhaps,  just  as  important,  whether  it  is  possible  to  tell 
when  not  enough  are  being  employed. 

Figure  1 is  a plot  of  the  product  of  the  rms  values  of  x and  PK  and  the  same  product  for  y 
and  Pv  of  a matched  K-V  distribution  in  a thin  lens  quadrupole  system  with  a 90°  phase 
advance.  The  space  charge  strength  is  sufficient  to  depress  the  tune  to  30”.  The  initial  condi- 
tions are  found  by  iterating  so  that  the  variation  in  x and  y of  the  bounding  ellipse  is  minimized 
over  the  first  pair  of  lenses.  This  seems  to  result  in  a favoring  of  the  x dimension  so  that  the 
rms  emittance  growth  in  the  x dimension  is  less  than  in  the  y direction.  It  is  important  to  note 
that  both  emittances  stop  evolving  after  about  40  lens  pairs  and,  except  for  some  fluctuations, 
are  essentially  constant. 


The  problem  of  collisions  in  particle  in  cell  codes  has  been  investigated  in  detail  in  con- 
nection with  plasma  simulations9.  In  this  run  16,000  particles  were  employed  on  a 128  x 128 
grid.  If  the  velocity  corresponding  to  the  maximum  x-velocity  of  the  initial  ellipse  and  the  ini- 
tial plasma  frequency  <nn  are  used,  the  parameter 
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can  be  defined  as  a characteristic  length  similar  in  sense  to  the  debye  length.  The  number  of 
particles  in  an  area  A2  would  then  provide  an  estimate  of  the  system  collision  time.  In  this  par- 
ticular simulation,  the  beam  area  is  approximately  130X2.  The  plasma  time  is  2n/<on  =3  doub- 
lets. The  collision  time,  rr,  is  therefore  approximately 


/V  __  N . . , . 

75o~50  d0“b"'s 


For  16000  particles  this  is  about  320  doublets. 


(2) 


A close  look  at  the  curves  in  Fig.  1 shows  that  both  the  x and  y rms  emittances  do  appear 
to  have  a very  slow  secular  growth  time  on  the  order  of  several  hundred  doublets.  By  estimat- 

2 U‘ 


ing  the  slope  for  this  curve  and  normalizing  to  the  initial  emittance 


|—  = 10  Vdoublets 
dT 


or  t ~ 103.  This  slow  secular  growth,  the  lack  of  coupling  between  the  xand  y rms  emittances 
and  the  approximate  calculated  value  of  collision  frequency  all  point  to  the  collisionless  nature 
of  this  simulation.  In  fact,  since  the  distribution  function  is  rapidly  changing  as  the  particles 
move  between  lenses,  and  the  final  distribution  is  far  from  a K-V  distribution,  only  an  order  of 
magnitude  estimate  of  the  collision  time  is  appropriate.  By  decreasing  the  number  of  particles 
it  should,  however,  become  possible  to  see  a change  in  behavior  when  the  run  time  becomes  of 
the  order  of  several  collision  times.  An  even  more  radical  change  should  occur  when  the  colli- 
sion time  becomes  comparable  to  the  growth  time  of  the  instability. 


Figure  2 is  a plot  of  the  same  simulation  run  with  1024  particles.  This  would  predict  a 
collision  time  of  about  2$  doublets.  Note  that  the  qualitative  behavior  of  the  emittance,  i.e.  the 
initial  constancy,  rapid  growth,  and  leveling  off,  is  the  same.  The  saturation  level  is  also  about 
the  same  as  in  Fig.  1 . However,  some  rounding  of  the  knee  of  the  curves  and  a greater  secular 
growth  in  the  emittance  are  apparent. 
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A further  decrease  of  a factor  of  two  to  S12  particles  has  a more  radical  effect.  The 
features  of  Figs.  I and  2 are  sharply  washed  out.  The  calculated  collision  time  for  this  run 
would  be  about  13  doublets. 

Figures  3 and  4 show  similar  curves  for  2S6  and  128  particles  respectively.  Little  evidence 
of  the  same  qualitative  behavior  as  in  the  collisionless  regime  is  still  evident. 

d(f/t„) 

Figure  5 is  a plot  on  log-log  paper  of  the  slope  vs.  N,  where  r is  the  period  of 

dr 

the  lens  doublet,  and  N is  the  number  of  particles.  Figure  4 has  a break  in  the  slope  at  about  t 
— 60  doublets,  so  both  slopes  are  plotted  for  N — 128  The  straight  line  is  a -45°  line  drawn 
assuming  the  rate  of  growth  in  emittance  is  proportional  to  the  inverse  of  the  number  of  parti- 
cles so  that  n — constant.  The  line  drawn  corresponds  to  a constant  of  5.  This  com- 

dr 

pares  to  a proportionality  of  about  50  obtained  for  the  ratio  /V/t,  in  Eq.  (2)  where  r,  is  the  col- 
lision time  found  previously.  It  should  be  noted  that  Figure  5 indicates  that  the  slope  of  the 
emittance  in  Figure  1,  which  is  the  run  with  16000  particles,  can  not  be  explained  by  numerical 
collisions. 

These  runs  have  dramatically  illustrated  the  effects  an  insufficient  number  of  particles  can 
have  on  the  behavior  of  the  beam  system  being  simulated.  Furthermore  it  appears  necessary  to 
have  a collision  time  much  greater  than  the  growth  time  of  the  instability.  In  the  absence  of  a 
well  defined  theory,  simulations  should  either  have  enough  particles  so  the  collisionless 
behavior  is  not  open  to  doubt  or  several  runs  must  be  made  in  order  to  test  the  numerics. 
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Fig.  3 - Time  evolution  of  the  same  system  at  Fig.  I with  512  particles.  Major 
differences  are  now  visible  in  the  behavior. 


Fig.  4 - Tima  evolution  of  the  system  with  256  particles  showing  increased  collislonal 
growth  rate  In  the  emittance  and  further  disruption  of  the  ordered  evolution  of  the 

regions  of  instability  and  saturation  fat  Fig.  i.  t ■_ 


d(«/«0)/dr 
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Fig.  5 - Time  evolution  of  the  jystem  with  1 28  ptxticlei  thowing 
even  itronger  collitional  effect*. 


